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1. INTRODUCTION 

The QCD Lagrangian contains seven parameters (apart from the vacuum angle 6) 
which are not determined by theoretical considerations: The gauge coupling constant g 
and the quark masses m^, one for each flavor (a = u,d,s,c,b,t). After renormalization, 
these constants become scale dependent: 

127r 

9 - = 33 _ (ln/iVA^)"Hl + • • •} 

ma ^ m,(/i2) = rha In^VA')""^!! + ...}• (1) 

(The dots stand for higher loop corrections.) Experimental study of these parameters is 
necessarily indirect, since quarks and gluons do not exist as free particles. The measure- 
ment of the running coupling constant as(/i^) in different processes and at different scales 
provides a nontrivial test of QCD at short distances. |l| Heavy quark (c,b,t) masses can 
- in principle - be investigated in the framework of heavy flavor quarkonia spectroscopy. 
Finally, as argued below, hadronic r decays seem to provide a unique source of experi- 
mental information on light (current) quark running masses mu,mfi,ms. The purpose 
of the present communication is to indicate how this information can be extracted in 
practice from experiment. 

Light quark masses (LQM) have a reputation of being "not well-measurable" 
and, indeed, their experimental determination has so far not even been attempted. On 
the other hand, there exists a huge number of theoretical estimates of their values, 

some of which claim an amazing accuracy. [^,[^ These estimates will be commented on 
shortly. At any rate, LQM are certainly tiny compared to the scale Ah ~ 1 GeV at 
which massive bound states of QCD are formed. This fact represents the major difficulty 
in the experimental study of LQM, since their contribution to hadron masses (except 
pions) is small and hard to estimate theoretically. On the other hand, LQM measure the 
absolute strength of chiral symmetry breaking. This follows from the fact that in QCD 
the divergences of observable axial and vector weak-transition currents are given by the 
equations 



d^idjtilbu) = [rrid + mu)di'j5U (2a) 

d'^is'jf.'j^u) = {nis + mu)si'j5U (2b) 

d^{d^^u) = {rrid - rriu) idu (2c) 

d^is-^^u) = {rris — niu) isu (2d) 



which are valid at all scales. The experimental determination of LQM is based on these 
equations, combined with the short distance properties of QCD. 

Before we develop the details of our argument that experimental investigation of 
LQM is feasible, let us briefly summarize why it is important: 

(i) It is important to know all parameters of the standard model. 
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(ii) The values of LQM are closely related [^] to the value of the quark- antiquark conden- 
sate (ipip), which is an important quantitative characteristic of the nonperturbative 
chiral order in the QCD vacuum. 

(iii) It is important to have a direct experimental check of the theoretical prediction for 
the ratios of LQM that is based on standard chiral perturbation theory. The 
existence of a possible disagreement between this prediction and some low-energy 
data has recently been pointed out. [p|,|lO| 



(iv) The claim that "m„ is not equal to zero" |TT[ should be framed in terms of a bound 
referred to a confidence level with a firm experimental basis. This issue is important 
for the understanding of the strong CP violation problem. 



2. HOW BIG IS THE DIVERGENCE OF THE AXIAL CURRENT? 

We shall mainly concentrate on the average of the u and d quark masses, 

m = {rriu + md)/2, (3) 



which, according to Eq. (2a), controls the strength of the divergence of the axial current. 
(The remaining cases of m,^ — and — will be briefly mentioned in Section 5.) 
The object of our concern will be the spectral function 

P(Q') = ^E(2^)'^^'HQ-^n)|(r^|9'^(rf7,75n)|0)|^ (4) 

271 „ 

where the sum extends over all states with the quantum numbers of the pion and with 
squared invariant mass Q^: 

n = vr", n^n'^n~ , tt^tt^tt^, n^n^n^n^n^ , .... (5) 



The spectral function p{Q'^) measures the amount of explicit chiral symmetry breaking 
at squared momentum transfer Q^. For large Q^, it is given by QCD perturbation 
theory: 

pm - i^.v^m? |i + f ^ +...}• (6) 

Hence, rh{Q'^) is directly measurable to the extent that p{Q'^) is measurable for sufficiently 
large Q^. The spectral function is comprised of individual exclusive components 

p{Q') = 2F'^Mt5{Q' - Ml) + p,^{Q') + p^rM) + pUQ') + ■■■■ (7) 



The one-pion contribution is the only one which is known; it is given by the pion mass, 
= 139 MeV, and by the pion decay constant, = 93.1 MeV. The remaining compo- 
nents, P3n,PKKny ^tc. are terra incognita of particle physics; they hide the experimental 
information on the quark mass m that we are looking for. 
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The only way to measure the unknown components of the spectral function Eq. (4) 
seems to be via exclusive hadronic r decays, such as r — > z/^ + Svr, Uj. + Svr, z/^ + ttKK, etc. 
(See Fig. 1.) In order to extract the desired information from these decays, one faces two 
distinct problems: 



The final hadronic state (say Sir) excited from the vacuum by the virtual W consists 
of J = 1 and J = waves. The J = 1 wave, which is large and resonant (cf. the 
Oi resonance), is a subject of current experimental study [Q; however, for our 
purposes, it is of little interest. P3n{Q^) is given by the square of the J = part of 
the corresponding amplitude, which is O(m^), i.e. too small to be directly measured. 
A model-independent solution of this problem is presented in the next section, p^^^ 
can be reconstructed from the interference of J=0 and J=l waves, which is 0(m), 
and can be measured even for rh as small as 7 MeV, provided one has large enough 
statistics. 

The second difficulty is related to the use of the asymptotic formula Eq. (6) in a 
not quite asymptotic region of Q^. Not only is limited by the r mass m,-, but 
the differential decay rate is strongly suppressed near = m?^. This difficulty may 
be resolved by using QCD sum rules [|I^ for the two-point correlator of the axial 
current divergence. In general, these sum rules can be put into the form 

/■oo 

m\s,) = H-\w, So) / dQ' w{Q\ sq) p{Q^). (8) 



Here, w denotes a positive weight function that selects contributions with Q"^ < sq. 
In practice, one uses w = exp{—Q^/so) or w = Q'^"'0{so — Q^). H{w,So) is then 
defined by the large behavior of the two-point correlator. It consists of a part 
given by QCD perturbation theory and of a nonperturbative part parametrized in 
terms of vacuum condensates. The theoretical uncertainty in H{w, sq) decreases 
with increasing sq. 



Sum rules of the type shown in Eq. (^ have been extensively used in the past to 
estimate quark masses. These studies show a reasonable stability with respect to different 
choices of the weight function w, and they are rather insensitive to the uncertainty in 
the nonperturbative part of H{w,so)- The main problem in these estimates is not the 
sum-rule technique itself but rather the complete absence of experimental information 
on the magnitude of p{Q^) beyond the one-pion contribution. Retaining just the pion 
contribution to the integral of Eq. (|), and using the positivity of individual components 
Eq. (7) of piQ"^), one finds a lower bound 

m(l GeV^) > (4 - 5) MeV. (9) 



We propose to exploit the sum rules Eq. (^ once more, and to improve the lower bound 
Eq. (9) by using the measured component psniQ"^) as input, r decays make it possible to 
work at higher sq {sq < m^), thus reducing the systematic error due to the nonperturbative 
part of H{w, sq). 
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3. RECONSTRUCTION OF pg^ FROM r-DECAY EXPERIMENTS 

There are two Svr contributions to piQ"^)'- ti^ti^ti^ and tt^tt^tt^. They are not 
related by any symmetry, and consequently they should be measured separately and 
added afterwards: 

P3AQ')=P-AQ')+Poo-iQ')- (10) 

The kinematic analysis of the corresponding r~ z/^Stt decays is identical in both cases; 
the difference resides merely in possible experimental difficulty in detecting neutral pions 
with sufficient efficiency. Even if the latter could not be attained, experimental infor- 
mation on p ^ would still be valuable. This is because isospin symmetry implies the 

inequality 

P3AQ')>Ip-AQ')^ (11) 

which could lead to a considerable improvement of the lower bound Eq. (9) even if no 
experimental information on were available. 

3a. Azimuthal asymmetries 

The kinematics of exclusive decays of r into three hadrons has been analyzed in full 
generality in Ref. |jT^ (hereafter referred to as KM). Here, we will use the same notation, 
except for a few simplifications which are appropriate for the special case of the decays 

T ^ 7r"(gi)7r"(g2)7r^(g3)i^., (12) 

where a = —(0), b = +(— ). These decays proceed via the axial weak current 

(7r"(gi)7r~(g2)7r+(g3)M7M75?^|0) = g2, ga) (13) 

and are described by three independent form factors. (The anomalous vector current 
contribution vanishes because of G-parity conservation.) Following KM, it is convenient 
to work in the hadronic center of mass frame, in which the three pions move in the x — y 
plane. The x axis is chosen parallel to qs and the hadronic plane is oriented so that the 
z axis points in the direction of qi x q2, where |gi| > \q2\. In this frame, the ^-component 
of Eq. (13) vanishes and the three independent components A^iAy and At are functions 
of the scalar hadronic variables 

q2 = + ^2 + ^3)2^ Si = (g - qif, S2 = {Q- q2f. (14) 

The three form factors A^, Ay and At are related to the form factors Fi, F2 and F4 of KM 
by 

A, = xiFi + X2F2, Ay = xs{F, - F2), At = Jq^F^, (15) 
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where the Xi arc kinematic functions defined in Eq. (33) of KM. The space components 
and Ay describe the J = 1 part of the matrix element. They are large and insensitive 
to the quark mass. The time component A^ is proportional to the matrix element of the 
divergence of the axial current, and consequently to the running quark mass rh. 

We do not assume a polarized r beam, and the knowledge of the direction of flight of 
the T in the lab system (i.e., of the r rest frame) is not needed. One should just measure 
the momenta of the three pions and reconstruct their center of mass frame as defined 
above. This determines Q^, the Dalitz plot variables si and S2, and the angles f3 (polar) 
and 7 (azimuthal) which define the direction of fiight Ul of the laboratory with respect 
to the oriented center of mass hadronic plane (see Fig. 2). The remaining variable is then 
9 (the r decay angle) given in terms of the hadronic energy and the energy E^. of the 
T, both in the laboratory frame. 

The hadronic structure functions of interest are 

Wab{Q\Ai)= [ dsids2Re{AlAi,), a,b = x,y,t, (16) 

JAi 

where Aj denotes a bin of the Dalitz plot. W^x^Wyy and W^y are independent of the 
quark mass, while the interference terms Wxt and Wyt are 0(m), and Wu is O(m^). The 
latter determines the contribution to the spectral function Eq. (4), 

P..M)-^,Y.Wu{Q^^i). (17) 



where the sum extends over all bins of the Dalitz plot. (By definition, the sum in (17) 
is independent of the manner in which the phase space is divided into bins.) The most 

straightforward measurement of the relevant structure functions Wab{Q'^ ^ ^i) involves the 
differential decay rate integrated over the polar angle /3 and over the r-decay angle 0: 



where 



3m^ 

W{Q\ A,) ^ WUQ\ A.) + Wyy{Q\ AO + ^—i—Wu{Q\ A,) (19) 



describes the differential decay rate integrated over all angles. The normahzed distribution 
in the azimuthal angle 7 is of the form 

/(7) = 1 + X2{A cos 27 + E sin 27) + Ai {Clr cos 7 + Cud sin 7) . (20) 



The coefficients A and B are "large," i.e., they are 0(1) in the chiral limit: 
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The small chiral symmetry breaking shows up through the left-right and up-down asym- 
metry coefficients 



_ TT 3ml W^t ^ TT 3m^ Wyt 

^^""^ 2m2 + 2g2 ly' 2ml + 2Q^W' ^ > 

which are proportional to the quark mass m. Measurement of these azimuthal asymme- 
tries represents the hard core of our method of determining rh. The coefficients Ai and 
A2 in Eq. (20) are kinematic functions of and of the velocity of the r in the labora- 
tory frame, — ^1 — rn^/El^ where Eq is the energy of the beam. The result from 
integration over the r decay angle 9, 

/I fl pos f) 
_^^-P„(cosV'), (23) 

where P„ are Legendrc polynomials and cosil) is a function of cos^, and (3^. as defined 
in KM. The shape of these functions depends strongly upon the energy of the machine, 
Eq. For r leptons produced at rest cos^ = 1, so one has 

A„(Q^0) = 1, (24a) 
whereas in the ultrarelativistic limit (3r — 1 these functions become (see Fig. 3.) 

Al(g^l) = 



{ml - Q'^y 

'ml — g^' 



X,{Q\ l) = -2 + 34^Al(g^ 1). (24b) 



For low-energy machines such as tau-charm factories, the integration over the decay 
angle 6 represents a simplification which does not lower the sensitivity to the azimuthal 
asymmetries (21) and (22). On the other hand, for high-energy machines such as at LEP 
{P^ = 0.999), CESR {/3r ~ = 0.93) or B-factories, Eq. (24b) indicates an important 
loss in sensitivity. In this case, one should make use of the knowledge of the decay angle 
distribution rather than simply integrating over it. 

3b. The trick 

One may take advantage of the smallness of the quark mass and neglect Wu com- 
pared to Wxx + Wyy in Eq. (19). Within this approximation, the measurement described 
above yields the structure functions Wxx, Wyy, Wxy and Wxt, Wyt for a given and for a 
given bin Aj of the Dalitz plot. On the other hand, similar information on the interference 
terms Im A'^At and Im A* At would require both a known nonzero r polarization and the 
reconstruction of the r rest frame. This information is, fortunately, not needed. We will 
prove that, provided the binning of the Dalitz plot is sufficiently fine, the O(m^) quantity 
Wtt can be reconstructed from the experimentally accessible 0(m) interference terms 
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0(g^A,)= f J^^M (25a) 



and from the spin-1 structure functions 



Given two complex numbers x and z, it is obviously impossible to reconstruct \z\^ 
from |a;p and Re{x*z) alone. However, for three complex numbers xi,X2 and z one has 
the identity 

2 

\z\^ = J2 Re{x*z){k-^)ijRe{z*Xj), (26) 

where denotes the inverse of the matrix kjj = Re{xiX*). This identity applies to the 
decay r —>■ StiUt-, which is characterized by three (complex) form factors. For each bin 
Aj, let us define the quantity 

Wu{Q',A,) = 0^(g^ A,)K-^(g^ A,)0(g^ A,) (27) 

given entirely in terms of observables (25a) and (25b). If a bin A shrinks to a point P of 
the Dalitz plot, one has 

Wa,{Q\ A) = Re{A:{P)A,{P)) e{A) (28) 

where 

e{A) ^ j^dsids2 ^ 0. (29) 

Consequently, in the small bin limit, the indentity (26) implies 

WuiQ', A) = WuiQ', A) + 0{e'). (30) 



It is now straightforward to reexpress this result in terms of the spectral function 
Eq. (17). For a given division of the Dalitz plot into N bins, define 

1 ^ 

p['Li{Q')^-^,j:'^u{Q^\), (31) 

i=l 

where Wu is defined by Eq. (27). Consider now the hmit N ^ oo such that e(A) ~ 
0{N~'^). The "true" spectral function is then obtained as 

p__,(g2) = hm pi^) (g^). (32) 
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In practice, however, one deals with a finite number of bins of a finite size, and it is 
important to analyze what can be concluded in this case. Independently of the size of A, 
Eq. (16) defines a scalar product which satisfies the Schwarz inequality 

A similar lower bound holds under the substitutions W^t and W^x ^Lxi where 

W^^ and W^^ are obtained by the replacement 

A^^A'^^ giA^ + g2Ay, (34) 

Qi, g2 being two arbitrary real constants. In this way, one generates a whole class of lower 
bounds, and one can then ask which one is the best. Remarkably, the answer leads to the 
reappearance of the expression Eq. (27); one easily finds 

r^^,^^^-WuiQ',A,). (35) 

(91,92) w^^ 

Consequently, for any finite size bin Aj, one has the lower bound 

Wtt{Q^A,)>Wtt{Q^A,), (36) 

which, according to Eq. (30), is saturated in the limit of small bins. For any splitting of 
the Dalitz plot into N bins, one can sum Eq. (36) over all bins and divide by 5127r^; by 

Eq. (31), the left hand side is simply p +(Q^), which is, of course, independent of binning, 

while the right hand side is precisely p^^^ (Q^) ■ One thus obtains the lower bound 

P-+(Q')>pi^U^')' (37) 

which is saturated with increasing number of bins. 

It is worth noting the independence of our method on the geometry chosen in the 
hadronic plane. This is a reflection of the fact that the expression Eq. (27) is invariant 
under 

<P^g<P, K gKg^, (38) 

where is a two-by-two real matrix, representing a general linear (nonsingular) real 
transformation of a two-dimensional space spanned by A^ and Ay. 

3c. Back to the sum rules 

Suppose that the above program has been completed for both r~ — > 7r~7r~7r''" inl- 
and —> -K^-K^-K^Vr decays. That is, the azimuthal asymmetries have been measured 
in bins of the Dalitz plot and the corresponding functions p _^ and Pqq_ have been 
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constructed according to Eq. (31) for < Nmax allowed by the available statistics. One 
can then return to the QCD sum rules (|]) and define the quantity 



In so/A2 
ln/iVA2 



24/29 



H-\w,so) 



dQMQ',So){2F^M'j{Q'-M:) 

+pL^_l(Q^)+pifJ(Q^)}. (39) 



One expects that niNifJ', ^o) will depend rather weakly on the choice of the weight function 
w, and on sq in a typical range 



2 GeV^ < So < ~ 3.18 GeV^, 



(40) 



especially for large N. For sq in the range (40), the running quark mass m(/i) obeys 

m(/x) > m7v(/i, So), (41) 

and this lower bound should saturate rather rapidly with increasing A^. The saturation 
as well as the variation of the right hand side of Eq. (39) with sq can be controlled 
experimentally. The latter variation can be considered as a source of systematic error 
arising from imperfections in the method. Another source of error, which is more difficult 
to estimate, is related to both perturbative (higher orders [15]) and nonperturbative 
(condensates, instantons |T^) uncertainties in the high-energy factor H{w,so). 

4. WHAT IS EXPECTED 

Light quark masses are the only entries in the Particle Data Group compilations 
|r^] that are not based on a measurement but on theoretical estimates. We first briefiy 
recall why a direct estimate of m is problematic without experimental information on 
P37r(Q^). Then we proceed to a model-dependent numerical study of the various steps 
of the method described in the preceding section. In particular, the resulting statistical 
error of a measurement of m(l GeV) will be estimated, using the maximum likelihood 
method. 

For definiteness, the finite energy version of QCD sum rules will be used here, closely 
following Ref. Eq. (||) then takes the form 



m^(so) 



47r2 
3^ 



l + R2iso) + 2C^{0,)/sl 



1 /■■^o 



(42) 



and similarly, the definition Eq. (39) becomes 



m^(/i. So) 



lllSo/A^\24/29 4^ 



l + R2{so) + 2C,{0,)/sl 



X 



(43) 



where the two-loop expression for -R2(so) as well as a discussion of the value of the 
dimension-4 condensate (74(04) can be found in Ref. 0. We will use the value for the 
QCD coupling constant as as determined in Ref. [Q. 
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4a. Chiral Perturbation theory 

Unlike m, tlie order of magnitude of the difference — m is fairly easy to estimate: 
From the size of S'?7y(3) symmetry breaking, one may infer 

(m, - m)(l GeV^) = 100 - 300 MeV, (44) 



in agreement with sum-rule results. [0 (A possible measurement of rus —m^ in r Ut-Ktt 
decay will be briefly mentioned in the following section.) The estimate Eq. ( ^4]) can be 
combined with the presumed value of the quark mass ratio r = ms/rh = 26 to conclude 
that m(l GeV^) = 4—12 MeV, in agreement with the lower bound (9). Actually, even this 
rather conservative and crude estimate of m is doubtful. The value r = 26 is based on the 
standard chiral perturbation theory which assumes that the quark antiquark condensate 
— {ipip) is large compared to F^rus- The latter hypothesis has no clear experimental or 



theoretical support, and in fact it need not be correct RlO] for the actual value of m 



The generalized chiral perturbation theory [10| which admits a lower value of —{ijjijj) does 



not fix the quark mass ratio r = ms/rh; any value 6.3 < r < 25.9 is consistent with the 
mass spectrum of pseudoscalar mesons, the lower bound for r arising from the condition 
of vacuum stability. Hence, accepting the estimate Eq. (H3) constrains m to the range 



4 MeV < m(l GeV^) < 50 MeV. (45) 



However, it is clear that finding m significantly higher than - say - 10 MeV would imply a 
considerably lower value oi —{ipip) and of r = ms/m than the standard chiral perturbation 
theory could support. (Independent, though indirect, measurements of the quark mass 
ratio r are possible in low-energy vr — vr scattering, _K'^4 decays pO[ and from observed 



corrections to the Goldberger-Treiman relations. 0) This illustrates once more why even 
an experimental lower hound on m would be of considerable interest. 

It is instructive to illustrate the variation of m and {ipip) within a simple model, 
in which the 37r contribution to the spectral function is described by a narrow = 0~ 
resonance - the vr': 

PsAQ') = 2F^'Mt, 5{Q' - Ml). (46) 



(Such a resonance does indeed exist for M^-/ ^ 1.3 GeV, but it is not narrow, jl^) The 
constant F-^/ - the tt' analogue of = 93.1 MeV - describes the coupling of the vr' to 
the axial current. It is proportional to m [F^^i = in the chiral limit) and it is expected 
to be small compared to F^. Its value is, however, unknown and there is no reliable 
model available to pin it down. (In particular, quark model estimates (e.g., Ref. [Q) 
of Ft^i are trustworthy only to the extent that the model would correctly describe the 
small breaking of chiral symmetry.) Considering F^-' as a free parameter, one may use the 
narrow vr' model Eq. (^) within the FESR Eq. (42) in order to investigate the sensitivity 
of m to Ft^i. Setting M^./ = 1.3 GeV and fixing sq from the higher moment sum rules (as 
explained in Ref. 0), one finds that as F^^i increases from 1.5 MeV to 15 MeV, m(l GeV) 
slowly rises from 7 MeV to 45 MeV. At the same time, the value of the qq condensate 
—F~^{tljilj) rapidly decreases from ~ 1.3 GeV to ~ 40 MeV. A priori, there is no reason 
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to exclude a value for Ft^i as large as 10 or 15 MeV. However, it is clear that, in this 
case, the tt' (or Stt) contribution to the quark mass would largely dominate over the single 
pion contribution. Moreover, the slow rise of m would not compensate for the drop of 
— {i/jiIj), so that —27h{ipip) would be considerably lower than F^M^. Such a violation of 
the Gell-Mann, Oakes, Renner formula is allowed and expected within the generalized 
chiral perturbation theory, [T^ and its experimental confirmation would be an argument 



in favor of the latter. 

It is conceivable that the tt' width cannot be neglected, and that a more realistic 
model of psniQ"^) is provided by a Breit-Wigner formula. However, even if the corre- 
sponding width were known, the absolute normalization of PsniQ"^) remains unspecified. 
Dominguez and de Rafael have proposed to normalize P3n{Q'^) by its low behavior 

1 A/f^ 

ps.m - (47) 



as given by chiral perturbation theory. In this way, they obtain rh{l GeV) = 7.8±1.0 
MeV, where the quoted error merely arises from the input data {MT^/,r.,ri) and from the 
uncertainty brought in by the sum rules. It does not involve the possible error in the 
absolute normalization of PsniQ"^) based on Eq. (^). 

Apart from doubts about using the ~ behavior of ps^r to fix its value at 
~ M|, ~ 1.69 GeV^, the main uncertainty in the above method of normalizing p^^^ 
resides in the formula Eq. (^) itself. This formula is obtained using the standard chiral 
perturbation theory in an experimentally unexplored domain of low-energy n — tt interac- 
tion, where results strongly depend on the value of the quark antiquark condensate (ipip). 
The generalized chiral perturbation theory, |T0[ which parametrizes this dependence in a 



model independent way, leads to a modification of Eq. (|47|): 



(48) 



Here, parameter introduced in ||I0[ which is a function of the quark mass ratio 

r = iris/rh, varying from 0^,^ = 1 (r = 25.9) to Otttt = 4 (r = 6.3), and which can be 
measured in low-energy vr — vr scattering. |10| Hence, PsniQ"^) cannot be normalized using 
its low-Q^ behavior, since the latter is only known up to a factor of 1-13.5. Pz-niQ"^) has 
to be determined experimentally. 



4b. A model-dependent theoretical experiment 

In order to show how the method above might work in practice, we have generated 
data using a model for the form factors of Eq. (15) and performed the analysis including 
an estimate of statistical errors. We do not expect any model to give better than an order 
of magnitude estimate of the observable quantities that we seek. The form factors Fi, F2 
are chosen to be exactly those of KM (cf. 0,|l^). For the J = form factor F4 we follow 
KM in assuming the dominance of the vr' resonance that decays into pvr; however, we use 
the minimal Tr'vrp coupling (which they do not), so that 
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(49) 



where = + 3M^ 




Ml, - g2 _ i^r,,(g2) 



^[Q2 - (M, + M^nQ^ - (M, - MJ2]/(4Q2). 



(50) 



The parameters M^/ and F^/ as well as the spin-1 function, Bp, are taken directly from KM 
without alteration. For simplicity, we take the pion mass to be zero. We have introduced 
a dimensionless parameter ^ which sets the scale of the J — form factor (the numerical 
factor of 35 is for convenience, with units GeV~^ understood here). The parameter ^ 
plays the same role here as the constant F^^' did in the narrow-resonance model described 
in the preceding subsection; that is, it is ^ which determines the contribution of the J = 
spectral function to the quark mass, and so it is an unknown quantity. We will see below 
that our normalization is chosen so that, within the framework of the model we have 
adopted, ^ is of order unity for values of m in the range (45). Note that the model treats 
the 7r~7r~7r+ and n^n^n^ modes similarly, so we may simply combine the effects of the 
two modes by absorbing all normalization into the single parameter ^. 

The first step in our experiment is to generate the data, i.e., compute the differential 
decay rate F((5^, Aj, 7) of Eq. (18). We do this for three cases, ^ = 0.5, 1 and 2. Using the 
form factors Fi, F2, F4 described above, we compute from Eq. (14) the functions A^, Ay, At; 
these are then used to determine from Eq. (16) the functions Wab, which then give F. For 
our three choices of ^, we show in Fig. 4 the asymmetry coefficients A, B,Clr,Cud of 
Eqs. (21), (22), where the Wab are integrated over the entire Dalitz plot. Note that the 
azimuthal asymmetry coefficients Ci^RjCud are on the order of a few percent or less, 
much smaller than A {B is accidentally small here, which is due to the choice of axes in 
the hadronic plane). This is expected, since Clr,Cud are 0{m), while A,B are 0(1) 
in the chiral limit. Similarly, we may compute the functions Wxx,Wyy,Wxy,Wxt,Wyt as 
functions of Q'^,si,S2. In Fig. 5 we show (for ^ = 1) the functions Wab integrated over the 
entire Dalitz plot; clearly Wu is negligible. 

Next, we study the lower bound of Eq. (43) for four choices of binning of the Dalitz 
plot: divide the square < Si, S2 < into 1, 4, 16, or 64 equal-sized squares; this gives 
a total oi N — 1, 3, 10 and 36 bins respectively for these four cases, since the region of the 
Dalitz plot is simply a triangular half of this square region. For a given choice of binning, 
we compute for each bin Aj the lower bound Wtt{Q^,^i) from Eq. (27). After summing 
over bins, we show (for ^ = 1) in Fig. 6 how J2iLi WuiQ"^, Aj) approaches WuiQ'^) as the 
number of bins grows. Then, we obtain p^^J^{Q^) from Eq. (31), and thus determine 
miv(A*, So) from Eq. (43). The latter is shown (for ^ = 1) in Fig. 7b as a function of Sq for 
various A^. (We take fi = 1 GcV.) Note that the bound is essentially independent of Sq 
for 2 GeV^ ^ sq ~ ~ 3.18 GcV^. The figure clearly indicates that the choice of ,^ = 1 
corresponds to a lower bound on m(/i) of about 14 MeV. Similar results are presented in 
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Fig. 7a for ^ = 0.5 (a lower bound of 7 MeV) and in Fig. 7c for ^ = 2 (a lower bound of 28 
MeV). We see that the bound is roughly proportional to ^ in this range of ^. For smaller 
^ the pion contribution to the spectral function eventually dominates, and the bound on 
m(/i) reduces to that of Ref. @, independent of ^. 

4c. Maximum likelihood estimate of error 

Apart from any systematic experimental error, there will be some error arising from 
finite statistics; this can be estimated by the maximum likelihood method. For N^rut 
events, the likelihood function is 

Nevt 

C{mN{lJ',So)) = ]J r(Xi;mAr(/i, So)), (51) 

i=l 

where Xj denotes the measured values of the phase space variables (5^,Si,S2,7 for the 
event i, and F is normalized to unity. For fixed /i = 1 GeV and fixed Sq, there is a one-to- 
one correspondence between the lower bound rhi^{fi, sq) and ^; therefore, the measurement 
of ^ is equivalent to the measurement of ruNifJ', ^o)- The best estimate for ^ is the value 
which maximizes the likelihood £, or equivalently ln£; the estimated standard deviation 
for the measurement of ^ is simply 

(52) 

Hence, the standard deviation for the measurement of the lower bound rh = TiiNifJ', so) is 

o'rh = (T^{dm/d^)] (53a) 

then, using Eq. (43) to compute {dm/dC,) [the ^ dependence of m is contained in p^^_((5^) 
and pi^l{Q% 

(Tra = (rh/ [1 - ml/m'^] , (53b) 

where rho is obtained from m by putting ^ = 0, i.e., including only the contribution of 
the pion to the integral of Eq. (43). For our model calculation, the standard deviation 
afh varies as l/\/W^t and is weakly dependent on sq and ^ (cr^ and m/S, are roughly 
independent of S,)- For example, from 250000 r STTUr events with ^ at the r- 
charm factory, one finds 0"^ ~ 1 MeV for 2 GeV^ < < wi^. The standard deviation am 
(in MeV) is displayed in the table below as a function of the number of bins and of C,'- 





N=l 


N=3 


N=10 


N=36 


N=oo 


^=0.5 


0.5 


0.6 


0.8 


0.9 


0.9 


^=1.0 


0.6 


0.7 


0.9 


0.9 


1.0 


e=2.o 


0.6 


0.7 


0.9 


1.0 


1.0 



evt 
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It is important to realize that am increases rapidly with /3r, i.e., with beam energy, 
because of the factor Ai (see Eq. (20) and Fig. 3). Thus, the design of an experiment 
to measure the quark mass must be optimized with respect to the competing effects of 
increasing r production and decreasing sensitivity as the beam energy is raised from 
threshold. 

In the preceding analysis, no use was made of any information coming from mea- 
surement of the decay angle 6. In the computation of am, before integration of 6, dT/d^ 
is proportional to cosip. If one makes no use of the measurement of 6, then one first 
integrates T over 9, so {dV /d^Y is proportional to (Ai)^ = (/ cos?/^)^. As pointed out in 
Sec. 3a, for = there is no loss in sensitivity from this procedure; however, this is no 
longer true for 7^ 0, and for ~ 1 (characteristic of LEP and CESR) we find that am 
is increased by about a factor of 4. So, for example, with 10^ r SnUr events with f3 ^1 
at CESR, am^Q MeV. 

On the other hand, if one makes use of the measurement of 6, then one integrates 
{dV / d^Y /T over 9 in order to estimate the sensitivity. One can see that this will lead to a 
much less dramatic degradation of sensitivity: ignoring the relatively weak 9 dependence 
of the 1/r factor, the 9 integral is proportional to / cos'^ ip = (2A2 + l)/3. A detailed 
study of this question is beyond the scope of the present work. 

5. OTHER COMBINATIONS OF LIGHT QUARK MASSES 

In principle, exclusive r decays allow the measurement of the divergences of all four 
currents d'^^'^^u, 57^75^, d'^^u, and s7^m which appear in Eqs. (2a-d). This suggests that 
the method described in detail in the previous section for the case + could be 
extended to other combinations of LQM: + m^, rrts — rriu and — m„. We shall 
briefly comment on each of these cases. 

5a. rUs + rUu 

Equation (2b) seems at first sight to allow a straightforward extension of the mea- 
surement of rrid + rriu to rris + rriu. The corresponding J = spectral function defined as 
in Eq. (4) receives contributions from single K state and the continuum starts with Rtttt. 
The corresponding component px-Kn of the spectral function could, in principle, be mea- 
sured in T —>■ Ut- + Ktih decays. The Cabibbo suppression of the latter may be partially 
compensated by a considerably larger ratio of signal to background (J = to J = 1), 
which is due to rris ^ rh. Unfortunately, this decay receives an anomaly contribution 
from the vector current in addition to the three usual axial- vector form factors A^., Ay and 
At. The presence of essentially complicates the reconstruction of Wu and of pKirn from 
observable quantities. Actually, this reconstruction turns out to be possible only using 
a polarized r beam. |]14[ (This example emphasizes once more the lucky circumstances 
which make the measurement of p3,r and of ma + rriu possible.) 

5b. rUs — rriu 

The difference rrts — rriu controls the divergence of the vector current s^^u - cf. 
Eq. (2d). There is no single-particle contribution to the corresponding J=0 spectral func- 
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tion and the continuum starts with the Kit state. The component Pktt{Q'^) can be mea- 
sured in exclusive r decays 

R-ji^ + Ur, K^n- + z/^ (54) 



described by two form factors representing analytic continuations of the well-known Kes 
form factors f±{Q'^)- The J = and J = 1 combinations of /± can be separated by 
measuring the polar angle between and the direction of the K (or vr) momentum in 
the hadronic center of mass frame. Notice that pxniQ^) ^ which is not expected to be 
as small as PsniQ"^), since ^ m - should now be measured directly. (Since there are 
only two form factors, it is impossible to reconstruct pKn from the s — p interference. 
The latter is, however, interesting in itself; it provides model-independent information on 
K — 71 phase shifts. [^) Finally, a lower bound for rris — mu is provided by the sum rules 
analogous to (8). 



5c. ma — rriu 

The leading contribution to the spectral function including the divergence (2c) 



of 



the vector current d'j^u comes from the rjn state. The kinematics of the decay 



n r] + Ur 



(55) 



is completely analogous to the decays (54). However, in this case, the J=0 form factor 
containing the information about — has no particular reason to be suppressed with 



respect to the J=l form factor. p2| The latter also vanishes as — — > 0, and is 
given by isospin breaking effects such as r/ — 7r° mixing. ||^ The decay (55) is obviously 
rare, but its observation and study would provide new information both about — niu 
and about the ratio (m^ — mu)/{ms — rh). 



6. CONCLUSION 

In the present paper, we have argued that it may be possible to obtain an ex- 
perimental determination of the running quark mass m = {rrid + itLu)/2 based on the 
measurement of angular asymmetries at the 1% level in the decay r — * z/^ + Stt. The hard 
core of our argument is the trick, described in Section 3b, which allows to reconstruct the 
dominant 3n component psniQ"^) to the order O(m^) spectral function p{Q^), associated 
with the divergence of the axial current, from the measurement of the order 0(m) angular 
asymmetries. 

Combining PsniQ"^) as extracted from experiment with QCD sum rules leads to a 
lower bound for the running quark mass rh{p) as shown by the inequality of Eq.(41). 
The appearance of only a lower bound arises from the following: We have restricted 
our attention to the exclusive contributions to p{Q^) coming from one and three pion 
intermediate states only. Other contributions, such as KKn or bir intermediate states, 
could also be studied experimentally, but they are expected to be small due to the lack 
of phase space. Also, in reconstructing psniQ"^) we worked with a finite number of bins in 
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the fixed Q"^ Dalitz plot, just as one does in practice. Again, the question of convergence 
with increasing number of bins can be addressed experimentally. 

In order both to study the different steps of the above method and to get an order 
of magnitude estimate of the quantities we are looking for, we have generated a set of 
data from a theoretical model described in Section 4. We wish to stress that the model 
on which these estimates are based should not be given any further significance beyond 
its illustrative purpose. We have found that a sample of about 250,000 background free 
r — > z/7- + 37r events is needed in order to reduce the statistical error on the lower bound for 
m down to 1 Mev. This error becomes then comparable to the theoretical uncertainties 
associated with the sum rules analysis. 

As it stands, our analysis is mainly suited for the low energy machines, such as 
a Tau/Charm Factory: In this case it is possible to integrate over the r decay angle d 
without loss of sensitivity for the measured angular asymmetries. If, on the other hand, 
one would like to use a high energy source, such as LEP, CESR, B-factories, the method 
described above would have to be adjusted and include the measurement of the r decay 
angle ^, in order to compensate for the loss of sensitivity. 

In the present study, we have focused on the mass combination m„ + m^, which is 
rather special, due to the particular role of G-parity. It also remains the most interesting 
combination from the point of view of its theoretical and phenomenological implications. 
In contrast, + m„ cannot be given a lower bound via the same analysis, because the 
vector current contributes an additional form factor to the corresponding r — > i^,- + Xtttt 
decays. On the other hand, the quark mass differences — m„ and m,^ — rriu can be 
investigated in the quasi-two-body decays r — > Ut+Ktt and in the rare decays r — > UT+rjir, 
respectively. 

We conclude that a high statistics sample of tagged, exclusive rdecays (r — > + 
7r~7r~7r"^, 7r°7r°7r~, KKn', ...) would qualify as an authentic light quark mass spectrometer. 
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FIGURES 



FIG. 1. Exclusive hadronic r decays. 



FIG. 2. The oriented center of mass hadronic plane. 
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FIG. 3. (a) The function Ai((5^,/3) as a function of Q"^ /m^ for selected values of (3. 



FIG. 3. (b) The function X2{Q^,P) as a function of for /? = 1. 



FIG. 4. (a) The asymmetry coefficients A,B as functions of (in GeV). 



FIG. 4. (b) The asymmetry coefficient Clr as a function of (in GeV), for ^ = 0.5, 1 and 2. 



FIG. 4. (c) The asymmetry coefficient Cud as a function of (in GeV), for ^ = 0.5, 1 and 2. 



FIG. 5. (a) Wxx-, Wxy, Wyy integrated over the entire DaUtz plot, for ^ = 1, as a function of 
g2 (in GeV). 



FIG. 5. (b) Wxt,Wyt,Wtt integrated over the entire Dalitz plot, for ^ = 1, as a function 
(in GeV). 



FIG. 6. How J2iLi WttiQ'^, Aj) approaches WuiQ'^) as the number of bins N grows. 



FIG. 7. (a) Lower bounds TfiNiji, sq) (in MeV) as a function of sq (in GeV^) and for selected 
values of A'^, for fixed M = 1 GeV and ^ = 0.5. 



FIG. 7. (b) Lower bounds mjv(/i, sq) (in MeV) as a function of sq (in GeV^) and for selected 
values of N, for fixed = 1 GeV and ^ = 1.0. 



FIG. 7. (c) Lower bounds rhi^{p, sq) (in MeV) as a function of sq (in GeV^) and for selected 
values of N , for fixed /J^ = 1 GeV and ^ = 2.0. 
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